Size of the emission source and collectivity in ultra-relativistic p-Pb collisions 
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Abstract 



The interferometric radii in the system formed in ultra-relativistic proton-lead collisions are investigated in a framework 
T— I -based on event-by-event 3+1 dimensional viscous hydrodynamics. We argue that the most central p-Pb collisions under- 
going collective expansion behave similarly to the peripheral nucleus-nucleus collisions. The interferometric observables 
can serve as signatures of the formation of an extended fireball. 
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[ The collective nature of the evolution of ultra-relativistic 
■nucleus-nucleus (A-A) collisions from Relativistic Heavy- 

ilon Collider (RHIC) to the Large Hadron Collider (LHC) 
energies has been well evidenced, in particular with 
such phenomena as the harmonic flow or the transverse- 
momentum dependence of the interferometric radii. The 
description of the intermediate evolution stage with rel- 

lativistic viscous hydrodynamics yields a successful quan- 
[titative prediction at the level of, say, 15% for the most 
•relevant observables and for a wide range of centralities c, 
including rather peripheral collisions up to c ~ 70% [lH3 • 
The application of this collective approach to the proton- 
nucleus, not to mention the proton-proton (p-p) collisions, 
■is more questionable {8] but very intriguing [9|, as some 
'features typical for collective phenomena have been ob- 
served in the highest-multiplicity proton-lead (p-Pb) and 
p-p collisions as well. 

■ The p-Pb collisions at \/s n n — 5.02 TeV have been 
investigated at the LHC [19n23j , with the original motiva- 
tion to study the initial state effects 2J| and, in particu- 
lar, the saturation in QCD. Important observables in that 
■respect are the two-particle correlations in relative pseudo- 
rapidity and azimuthal angle. The long-range correlations 
in rapidity are formed in the very early stage of the colli- 
sion and present a dedicated probe of the initial state, on 
the other hand, the azimuthal correlations can be signif- 
icantly modified by the final state rescattering. When a 
fireball of strongly interacting matter is formed, azimuthal 
correlations due to collective flow appear in the interaction 
region 2^ 2^. For instance, in Ref. 27 1 we have argued 
that the appearance of the same-side ridge structure in 
the correlation data in the p-Pb collisions measured at the 
LHC is semi-quantitatively described with event- by- event 



hydrodynamics, which provides a strong case for the inter- 
pretation based on collective harmonic flow. An alterna- 
tive, appealing explanation of the appearance of the same- 
side ridge in p-Pb interactions is based on the color-glass 
condensate approach for the initial state 28, 2^. There- 



fore, it is important to look for independent estimates of 
the role of flnal state interactions in the dynamics. 

In this Letter we argue that the behavior of the pio- 
nic interferometric radii in the most central p-Pb colli- 
sions should be used as a flngerprint of collectivity; if the 
experimental results follow the pattern of A-A collisions, 
then collective behavior is present. This concerns both 
the values of the radii as well as their dependence on the 
transverse momentum of the pair, governed by the flow 
(30| . We evaluate the pionic interferometric radii for the 
most central p-Pb system with relativistic hydrodynam- 
ics and predict that in this treatment the system should 
closely imitate the peripheral nucleus-nucleus collision. If 
this is confirmed experimentally, it should serve as another 
strong case for the presence of collectivity in the most cen- 
tral p-Pb collisions. At the same time, the low-multiplicity 
p-Pb collisions are in our view not expected to display the 
above advocated behavior and should follow the p-p pat- 
tern, hinting different dynamics. Typically, the size and 
life-time of the collective source formed in central p-Pb 
collisions is 3-4 fm 
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The initial geometry of a small-source system formed 
in p-Pb collisions is dominated by fiuctuations, therefore 
the costly machinery of event-by-event viscous hydrody- 
namics must necessarily be applied to properly describe 
the azimuthally asymmetric components of the collective 
flow [32I. As the result, the hydrodynamic expansion in 
most central p-Pb collisions generates a sizable elliptic and 
triangular flow 31| , while the two-dimensional correlation 
functions in relative pseudorapidity and azimuth are in 
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Figure 1: Pion interferometry radii -Rsidc (e')i ^?out (b), and -Riong (c) for average pair momentum fcy = 400 MeV for different collisions systems 
and energies, plotted as functions of the charged particle multiplicity. The compilation of the experimental data is taken from 10, 11] . The 
data for the RHIC energies are from the STAR Collaboration [mUM and for the LHC energies from the ALICE Collaboration 14, 15] . The 
lines are to guide the eye. Various hydrodynamic calculations come from [ifiMlSH . The results of this work for the p-Pb system are indicated 
with filled stars (standard source) and encircled stars (compact source). 
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Figure 2: The model predictions for the pionic interferometric radii 
in p-Pb collisions for centrality 0-3.4%, plotted as functions of the 
average transverse momentum of the pair for the standard source 
(solid lines) and compact source (dashed lines). 



semi-quantitative agreement with experimental observa- 
tions ,21.-23.1 ■ following the same mechanism as in A- A 
collisions |33l. |34 1 . 

Interferometric correlations explored in this work are a 
snapshot of the emission points at the final stage 3^ 36 1. 
Fitting the Bertsch-Pratt [H, 111] formula to the same-sign 
two-pion correlation functions yields estimates for the fem- 
toscopic size of the emission source. Quite remarkably, the 
systematics of the interferometry radii in nuclear collisions 
at different energies shows an approximate extended scal- 
ing with the multiplicity of the system 11 , 1^ , 
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as shown in Fig. [T] We note that for the A-A collisions 
the hydrodynamic results shown in Fig. [T] are consistent 
with the data, properly reproducing the scaling ^ found 
experimentally. This stems from a generic relation be- 
tween the interferometry radii and the size of fireball in 
A-A collisions and is found in many hydrodynamic calcu- 
lations [39-42]. Admittedly, more precise probes, sensitive 
to the details of the collective flow profile, such as the ra- 
tio i?sido/-Rout, make the exact agreement more difficult to 
achieve in hydrodynamics, however, with a proper choice 
of the initial condition and equation of state [391441 1 it can 



also be accomplished. The cascade model approaches also 
work properly for the A-A system 43, 44]. 

On the other hand, the interferometry radii measured 
in p-p collisions show a similar scaling trend but with a 



very different slope than the A-A case (cf. Fig. [T]). This 
indicates that the mechanism responsible for the forma- 
tion of the interferometric correlations may be distinct in 
elementary 45| and nuclear collisions. We note that the 
hydrodynamic modeling of the p-p system assuming col- 
lective expansion and rescattering in the final state is not 
always compatible with the data 3 46-4^. We should 
admit that the analysis meets some difficulties: the ex- 
tracted interferometric radii in the small p-p system de- 
pend strongly on the form of the fitted correlation func- 
tion and the background subtractions, they require preser- 
vation of the conservation laws, moreover, depend on the 
resonance contributions or the effects of the uncertainty 
principle 14, \^ 47, [H^]. Another important ingredient 



in small systems that may influence the interferometry 
radii and their momentum dependence is the unknown pre- 
thermal flow [4lJ,|47[. All in all, the hydrodynamic predic- 
tions for the p-p pionic interferometric radii can overshoot 
the data by more than 50%, showing that the hydrody- 
namic description does not describe the p-p data in a uni- 
form way. 

The situation is hopefully different for the p-Pb system, 
where observed two-particle correlations [2l|-|23| suggest 
the possible existence of collective flow. If this picture 
is true, it means that the system formed in p-Pb colli- 
sions is sufficiently large and long-lived to accommodate 
a hydrodynamic expansion stage. In our model, the ini- 
tial condition is generated with GLISSANDO [51. The 
parameters of the calculations are similar as in [31J, ex- 
cept that they are adjusted for the collisions energy of 
y^SNN = 5.02 TeV. Thus we take the NN inelastic cross 
section ctnjsi = 67.7 mb, moreover, we use a realistic 
(Gaussian) wounding profile [H?] for the NN collisions. 
When the individual NN collision occurs, a source is pro- 
duced, meaning deposition of energy and entropy in a lo- 
cation in the transverse plane and spatial rapidity. In the 
conventional wounded nucleon model it is assumed that 
the sources are located in the transverse plane in the cen- 
ters of the participating nucleons. This amounts to rather 
large initial transverse sizes in the p-Pb system. Locating 
the source in the center-of-mass of the NN systerrQ is also 
admissible, which leads to a more compact initial distribu- 
tion. We use both variants, labeled standard and compact, 
which allows us to provide upper and lower limits for the 
size of the initial source and thus estimate the model un- 
certainty. The average rms radii for the two sets of initial 
conditions are 1.5 and 0.9 fm, respectively. The proba- 
bilistic nature of the Glauber model leads to initial source 
distributions which fluctuate event-by-event, referred to 
as the geometric fluctuations. No other possible sources 
of fluctuations in the initial phase, such as fluctuations of 
the color fields at smaller scales, are incorporated. 

In our simulations we use the event-by-event 3 + 1- 
dimensional viscous hydrodynamic model as described in 
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Figure 3: The model predictions for the central sections, in direction 
(Jsidci Qout and giong (panels a), b) and c) respectively), of the inter- 
ferometric correlation functions for the pion pairs (points) and the 
corresponding three-dimensional Gaussian fits (lines) for the most 
central (c = — 3.4%) p-Pb collisions. 



detail in 271, |3l[. The statistical emission and reso- 
nance decays at freeze-out are carried out with THER- 
MINAT0R2 The centrality 0-3.4% for the highest- 

multiplicity p-Pb collisions is defined in a simplified way 
by the condition on the number of participant nucleons, 
Npart > 18. We generate 450 distinct hydrodynamic evo- 
lutions of different initial conditions, both for the standard 
and compact source scenario. For each freeze-out hyper- 
surface at Tf = 150 MeV we generate 1000 THERMINA- 
TOR events to increase the statistics. The expansion of 
the compact source is faster than for the standard source. 
The stronger flow in the compact source case causes an 
increase of the average transverse momentum by 20%, but 
the size of the fireball at freeze-out is similar for the two 
considered scenarios. 

The femtoscopic correlation functions for the pion pair 
are obtained with the two-particle method described in 
detail in Ref. [Hi] and implemented in THERMINA- 
T0R2 Technically, the correlation functions involve 
pairs from the same hydro events in the numerator, and 
mixed pairs from different hydro events in the denominator 
of the correlation function 



Ciq,k) 



J d^xid''x2S{xi,pi)Six2,P2)\'^{k*,r*)f^ 



(2) 



/ d^XiS{Xi,Pi) J d'^X2S{x2,P2) 

where q — pi — p2 is the relative momentum of the pi- 
ous, k = {pi +p2)/2 is the average pair momentum, and 
the asterisk indicates the variables boosted to the pair rest 
frame [Hi]. For all the THERMINATOR events generated 
from the same hydrodynamic event the emission source 
S{x,p) is the same, hence pairs from such 1000 events can 
be combined in the model calculation to improve statistics 
in the numerator. Still, by combining pairs from different 
freeze-out hypersurfaces, or using the event- by- event aver- 
aged initial conditions, we find that the effects of the event- 
by-event fiuctuations of the emission source are small 17[ . 
The interferometric radii are obtained by fitting the Gaus- 
sian shape to the correlation functions. The Coulomb ef- 
fects, expected to be very small in the p-Pb system, are 
not taken into account in the pair wave function ^E*, and 
consistently, no Coulomb corrections are used in the fitting 
procedure. 

The result displayed in Fig. [T] shows that our hydro pre- 
dictions for the most central p-Pb system fall close to the 
A- A line, thus displaying the collective behavior. The dis- 
tinction from the p-p trend is clearly seen, in particular 
for the standard-source case for i?out and i?iong, where the 
difference is about a factor of 2. The compact source leads 
to somewhat smaller interferometric radii, in particular for 
^?out which is reduced by 25%. The other radii are very 
little affected by the initial source size. 

In Fig. [5] we give the dependence of the pionic femto- 
scopic radii on the transverse momentum of the pair, kx- 
Again, we show the case of the most central p-Pb colli- 
sions, as these are most likely to display collectivity. We 
note the expected fall-off of the radii with kr caused by 
the collective flow. For the compact initial distribution the 
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values of the radii are above 1.8 fm for i?sido and above 2 fm 
for i?out and i?iong- For the compact case the i?out radius 
is visibly reduced, while i?out and i?iong are only slightly 
smaller. Observing experimentally such large sizes (com- 
pared to the proton radius) in the most central p-Pb would 
demonstrate the formation of a fireball. 

In femtoscopic studies the shape of the correlation func- 
tion and its departure from the Gaussian form is frequently 
studied 47|, as it carries relevant information on the dy- 
namics of the system. We provide this information for 
the standard source in Fig. |31 where we plot the central 
sections of the correlation function. The departure from 
the three-dimensional Gaussian fit is clearly visible. In 
particular, at low momenta qa the correlation function is 
much sharper, rising well above the fitted Gaussian profile. 
Nevertheless, following the experimental works reported in 
Fig. [T]we use the three-dimensional Gaussian profile in ex- 
tracting the interferometry radii. 

To conclude, we state again the importance of the ex- 
perimental femtoscopic measurements for the p-A system, 
which will help to determine the nature of its dynamics. 
The proximity to the A-A scaling line of Fig. [1] will place 
the system in the collective evolution mode, on the other 
hand, if it turns out to be close to the p-p line, elementary 
dynamics will be vivid. Our simple hydrodynamic calcula- 
tion for the most central p-Pb system gives radii consistent 
with the A-A scaling. We should note, however, that in 
a more realistic treatment we expect some deviations due 
to remnants from the elementary p-p collisions, as mod- 
eled for instance in the core-corona picture. We also note 
that if a large size fireball is found, it could be used in 
quenching models to be compared with the Raa data. 
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